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Abstract 

We consider the Rayleigh- Taylor problem for two compressible, immiscible, inviscid, 
barotropic fluids evolving with a free interface in the presence of a uniform gravita- 
tional field. After constructing Rayleigh- Taylor steady-state solutions with a denser 
fluid lying above the free interface with the second fluid, we turn to an analysis of the 
equations obtained from linearizing around such a steady state. By a natural varia- 
tional approach, we construct normal mode solutions that grow exponentially in time 
with rate like e*v / l£l, where £ is the spatial frequency of the normal mode. A Fourier 
synthesis of these normal mode solutions allows us to construct solutions that grow 
arbitrarily quickly in the Sobolev space H k , which leads to an ill-posedness result for 
the linearized problem. Using these pathological solutions, we then demonstrate ill- 
posedness for the original non-linear problem in an appropriate sense. More precisely, 
we use a contradiction argument to show that the non-linear problem does not admit 
reasonable estimates of solutions for small time in terms of the initial data. 

1 Introduction 

1.1 Formulation in Eulerian coordinates 

This paper concerns the compressible, inviscid Rayleigh- Taylor problem in an infinite slab 
Q := M 2 x (—m, t) C M 3 . For this we consider two distinct, immiscible, inviscid, compressible, 
barotropic fluids evolving within Q for time t > 0. The fluids are separated from one 
another by a moving free boundary surface E(t) that extends to infinity in every horizontal 
direction; this surface divides fl into two time-dependent, disjoint, open subsets fl±(t) so 
that D, = Q + (t ) U (t ) U E (t) and E (t ) = fi+ (t) H (t) . The fluid occupying fi+ (t) is called 
the "upper fluid," and the second fluid, which occupies fi-(£), is called the "lower fluid." 
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t Supported by an NSF Postdoctoral Research Fellowship 



1 



The two fluids are described by their density and velocity functions, which are given for each 
t > by 

p±(-,t) : Sl±(t) ->- R + = (0,oo) and u±(;t) : Q±(t) ->■ ffi 3 (1.1) 

respectively. We shall assume that at a given time t > the density and velocity functions 
have well-defined traces onto £(i). 

We require that the fluids be sufficiently smooth to satisfy the pair of compressible Euler 
equations: 

d t p± + div(p±u±) = for t > 0,x G fi±(t) 

p±(d t u± + u±-Vu±)+V(P±(p±)) = -gp±e 3 for i > 0, x G Sl±(t). 

Here we have written g > for the gravitational constant, e3 = (0, 0, 1) for the vertical unit 
vector, and — ge^ for the acceleration due to gravity. We have assumed a general barotropic 
pressure law of the form P± = P±{p) > with P± G C°°((0, oo)) and strictly increasing. We 
will also assume that 1/P±. G L^ c ((0, oo)). Finally, in order to produce the Rayleigh- Taylor 
instability, i.e. construct a steady state solution with an upper fluid of greater density at 
£(£), we will assume that 

Z := {z G (0, oo) | P_(z) > P + (z) and P_(z) G P+((0, oo))} ^ 0. (1.3) 

In particular this requires the pressure laws to be distinct, i.e. P_ ^ P + . For a general 
discussion of the physics related to these equations and to the Rayleigh- Taylor instability, 
we refer to [5] and the references therein. 

A standard assumption is that both the normal component of the velocity and pressure 
must be continuous across a free boundary between two inviscid, immiscible fluids [Sj. This 
requires us to enforce the jump conditions 

■ u + )\ m ~ ■ u„)\ m = 

(p + (p+))M-(P-(p-))M = o, 

where we have written the normal vector to £(i) as v and /|s(t) for the trace of a quantity / 
on S(t). We will also enforce the condition that the normal component of the fluid velocity 
vanishes at the fixed upper and lower boundaries; we implement this via the boundary 
condition 

U-(xi,X2, —m, t) ■ 63 = u + (xi, £2, t) ■ e 3 = for all (xi, X2) G M 2 , t > 0. (1.5) 

The motion of the free interface is coupled to the evolution equations for the fluids (11.21) 
by requiring that the surface be advected with the fluids. This means that the velocity of the 
surface is given by (u ■ v)v. Since the normal component of the velocity is continuous across 
the surface there is no ambiguity in writing u ■ v. The tangential components of u± need 
not be continuous across £(£), and indeed there may be jumps in these. This allows for the 
possibility of slipping: the upper and lower fluids moving in different directions tangent to 
£(t). Since only the normal component of the velocity vanishes at the fixed upper and lower 
boundaries, {x 3 = £} and {£3 = — m}, the fluids may also slip along the fixed boundaries. 
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To complete the statement of the problem, we must specify initial conditions. We give 
the initial interface X , which yields the open sets f2±(0) on which we specify the initial data 
for the density and velocity, p±(0) : f2±(0) — > 1R + and u±(0) : f2±(0) — > M 3 , respectively. 

It will be convenient in our subsequent analysis to rewrite the second Euler equation 
using the enthalpy function 

r z P' (r) 

h±(z) = / -^^dr. (1.6) 
J\ r 

The properties of P± guarantee that h± G C°°((0, oo)) are both strictly increasing, and hence 
invertible on their images. Then the second compressible Euler equation may be rewritten 

as 

d t u± + u± ■ Vu± + V(/i±(p±)) = -ge 3 . (1.7) 



1.2 Steady-state solution 

We seek a steady-state solution with u± — and the interface given by {x 3 = 0} for all 
t > 0. Then fi+ = Q + (t) = R 2 x (0,£),Q- = Sl-(t) = R 2 x (-m,0) for all t > 0, and the 
equations reduce to the ODE 

d(h±(p±)) 

= -g m fi ± (1.8) 

dx 3 

subject to the jump condition 

P + (p + ) = P_(P-)on{x3 = 0}. (1.9) 

Such a solution depends only on x 3 , so we may consolidate notation by assuming that p± 
are the restrictions to (0,£) and (— to, 0) of a single function p = po(x 3 ) that is smooth on 
(—m, 0) and (0, £) with a jump discontinuity across x 3 = 0. 

To solve this we recall that h± e C°°((0,oo)) are both strictly increasing and invertible 
on their images. The solution to the ODE is then given by 

fhZ 1 (h-(po)-gx 3 ), -m<x 3 <0 
Po{x 3 ) = < ^ + (1.10) 

(/i+(Po ) - ffx 3 ), o < x 3 < e, 

where p$ > is a free parameter satisfying P_(p ) e P + ((0, oo)), which allows the jump 
condition to be satisfied by choosing Pq > according to 

p+ = P^(P-(Po ))• (1.H) 
For p to be well-defined, we will henceforth assume that £, m > are chosen so that 

(MPo ) + ^) e /i-((0, oo)) and (/i + (p+) - ^) e ^+((0, oo)). (1.12) 

Note that p is bounded above and below by positive constants on (—m,£) and that p is 
smooth when restricted to (— m, 0) or (0, £). 

Since we are interested in Rayleigh- Taylor instability, we want the upper fluid to be 
denser than the lower fluid at the interface, i.e. Pq > Pq . This requires us to choose Pq so 
that 

P+\P-(Po)) > Po P-(Po ) > P + (Po)- (1-13) 
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The latter condition is satisfied for any p £ Z, where Z was defined by (II. 3p ; we assume 
Pq takes any such value. Then 

[p ]:=p+-Po>0. (1.14) 

For the sake of clarity, we include an example of the solution, po, when the pressure laws 
correspond to polytropic gas laws, i.e. P±(p) = K±p y± for K± > 0, y± > 1. The solution is 
then given by 



Po{x 3j 



(Po 



_-i _ g(7--i) 

if_7_ 



1/(7—1) 



x 3 < 



(Po + 



^- V/(7+_1) 0<x 3 <^ TT (p + V-- ! 







9( 7+ -l)^0 
X 3 ^ 9 ( 7+ -l) [PO ) 



;i.i5) 



with modification to solutions po(xs) = p$ exp(— gxs/K±) when either 7+ or 7_ is 1. The 
jump condition requires that 



K. 



1/7+ 



K-) (Po 



\ 7-/74 



(1.16) 



For a polytropic gas law, the condition that p$ > p is equivalent to 

K \ K 



1.17) 



If 7+ = 7_ this requires K_ > K + and any choice of p > 0. If 7+ 7^ 7_ then K + > 
can be arbitrary, but we must require that Pq > satisfies 



A' 4 



Po < (i?7 



l/(7-~7+) 



1/(7+ -7-) 



if 7_ > 7 4 



1.18) 



if 7+ > 7- 

In either case, to avoid the vanishing of po, £ is chosen so that 



< i < - 7 ^ ± T t(p T + -\ 

9{1+ - 1) 

but the parameter m > may be chosen arbitrarily. 



;i.i9) 



1.3 Horizontal Fourier transform and piecewise Sobolev spaces 

Before stating the main results, we define some terms that will be used throughout the paper. 
For a function / G L 2 (Q), we define the horizontal Fourier transform via 
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By the Fubini and Parseval theorems, we have that 



\f(x)\ 2 dx 



An 2 



d^dxz- 



1.20) 



1.21) 



We now define a function space suitable for our analysis of two disjoint fluids. For a 
function / defined on Q we write / + for the restriction to Q + = M? x (0, £) and /_ for the 
restriction to = IR 2 x (— m, 0). For s£i, define the piecewise Sobolev space of order s 
by 

H'(Q) = {f\f+e H S (Q + ), /_ G H S (Q_)} (1.22) 

endowed with the norm = + For ^^ ffe can take the norms 

to be given as 




d£ (1.23) 



for I = (— m, 0) and i+ = (0, £). The main difference between the piecewise Sobolev space 
H S (Q) and the usual Sobolev space is that we do not require functions in the piecewise 
Sobolev space to have weak derivatives across the set {£3 = 0}. 



1.4 Summary of main results and plan of paper 

The motion of the free surface £(t) and the domains fl±(t) present several mathematical 
difficulties, so we begin our analysis by switching to a Lagrangian coordinate system in 
which the domains of the upper and lower fluids stay fixed in time as Q + = IR 2 x (0, t) 
and Q- = R 2 x (— m, 0), respectively. Since the steady-state solution has vanishing fluid 
velocity, the steady-state constructed above in Eulerian coordinates is also a steady-state in 
Lagrangian coordinates. 

The first part of the paper is devoted to a study of the equations obtained by linearizing 
the compressible Euler equations, written in Lagrangian coordinates, around the steady- 
state solution. The resulting linear equations have coefficient functions that depend only on 
the vertical variable, X3 G (—m,£). This allows us to seek "normal mode" solutions (cf. [TJ) 
by taking the horizontal Fourier transform of the equations and assuming the solution grows 
exponentially in time by the factor e A ®*, where £ G IR 2 is the horizontal spatial frequency 
and A(£) > 0. A similar strategy was employed in the non-barotropic, horizontally periodic 
case in [7]. This reduces the equations to a system of ODEs for each £, which constitute a 
variant of a classical boundary- value eigenvalue problem with eigenvalue A(£) since one of the 
unknown functions in the system only appears with first order derivatives in the equations. 

In spite of its non-standard structure, the eigenvalue problem is amenable to solution by 
constrained minimization, which we then employ. The use of variational methods to produce 
the normal mode solution for each £ is essential to our analysis since it gives rise to detailed 
estimates of the behavior of A(£) as £ varies in IR 2 . We show in Lemma [3.51 that there is a 
constant C so that 

0<C<^< v ^as \C\^oo. (1.24) 
V Isl 

Since A(£) — > 00, normal modes with a higher spatial frequency grow faster in time, which 
provides a mechanism for Rayleigh- Taylor instability. Indeed, we can form a Fourier synthesis 
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of the normal mode solutions constructed for each spatial frequency £ to construct solutions 
to the linearized compressible Euler equations that grow arbitrarily quickly in time, when 
measured in H k (Q) for any k > 0. This is the content of Theorem 13.101 The variational 
methods used are quite general and robust; indeed, in [3] we extend the methods of this 
paper to construct growing solutions to the linearized Navier-Stokes equations with surface 
tension and density-dependent viscosity coefficients. 

As a preliminary for our analysis of the well-posedness of the non-linear problem, we then 
consider the well-posedness of the linearized problem. Inspired by a result in [3] , we show 
a connection between the growth rate of arbitrary solutions to the linearized equations and 
the eigenvalues A(£), which then gives rise to a uniqueness result, Theorem 14.41 In spite of 
the uniqueness, the linear problem is ill-posed in the sense of Hadamard in H k (Q) for any 
k since solutions do not depend continuously on the initial data. This is shown in Theorem 
14.51 by employing Theorem 13.101 to build a sequence of solutions with initial data tending 
to in H (ft), but which grow to be arbitrarily large in H k (Q) arbitrarily quickly. Again, 
the construction depends heavily on the detailed knowledge of the normal mode solutions 
provided by the variational methods. 

With linear ill-posedness established, we then prove ill-posedness of the fully non-linear 
problem. There is no general theory that guarantees the ill-posedness of a non-linear problem 
given the ill-posedness of the resulting linearized problem. As such, it is novel and remarkable 
that in the present case, linear ill-posedness does indeed give rise to ill-posedness for the 
non-linear problem. This is in some sense a compressible analogue to the Rayleigh- Taylor 
ill-posedness results in the incompressible regime, proved in [2]. 

To see ill-posedness, we rewrite the equations in Lagrangian coordinates with the un- 
known functions given as perturbations of the steady-state solution (see ( 15. 5 p for the exact 
formulation). For any k > 3, we then define a notion of non-linear well-posedness, which we 
call property EE(k) (see Definition 15.11 for the precise statement). Property EE(k) requires 
local-in-time (on an interval (0,t )) existence of solutions for initial data with small H k (Q) 
norm, along with L°°((0, to); H 3 (Q)) estimates of the solutions (written here generically as 
X(t)) of the form 

sup \\X(t)\\ H3(n) <F(\\X(0)\\ Hk{n) ) (1.25) 

for some function F, satisfying the Lipschitz condition F(z) < Cz for some C > for all z 
in a neighborhood of 0. Condition EE(k) is quite general and a reasonable choice for any 
well-posedness theory. 

We then show in Theorem 15.21 that it is impossible for property EE(k) to hold for any 
k > 3, which implies ill-posedness for the non-linear compressible Euler equations. The 
important feature of EE(k) is that k > 3 is arbitrary. If the initial data are extremely 
smooth (k very large), the failure of property EE(k) means it is impossible to control even 
the H 3 (Q) norm for small time. The proof of Theorem 15.21 is a reductio ad absurdum. 
Indeed, we show that if EE(k) holds, then it is possible to obtain certain estimates for the 
corresponding linearized equations that violate our linear ill-posedness result, Theorem 14.51 

All of our analysis is performed on the semi-infinite slab Q = R 2 x (— m, £), and one 
may wonder if the Rayleigh- Taylor instability is somehow related to the infinite extent of 
the horizontal component R 2 . This is not the case; all of our results carry over essentially 
word for word if Q is replaced by Cl = T 2 x (— m,£), where T 2 = S 1 x S 1 is the 2— torus. In 
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this case, all of the functions are required to be periodic in the horizontal directions, and 
the horizontal Fourier transform (11.201) with continuous spatial frequencies £ G M? must be 
replaced with the horizontal Fourier transform of T 2 , for which the spatial frequencies are 
constrained to £ e Z 2 . For the sake of brevity we will not rewrite the results for Q. There 
is also nothing essential in our analysis about the domain being contained in IR 3 . The same 
results hold if Q or Q is replaced by Q = M 1 x (— m, £) or Q = T 1 x (— m, £). 

The paper is organized as follows. In Section 2 we perform the switch to Lagrangian 
coordinates and record the linearized equations. In Section 3 we construct the growing 
solutions to the linearized equations. In Section 4 we analyze the linear problem, proving 
uniqueness and discontinuous dependence on the initial data. In Section 5 we prove the 
ill-posedness result for the non-linear problem. 

2 Formulation in Lagrangian coordinates 
2.1 Lagrangian coordinates 

The movement of the free boundary and the subsequent change of the domains fl±(t) in 
Eulerian coordinates create numerous mathematical difficulties. We circumvent these diffi- 
culties by switching to Lagrangian coordinates so that the free interface and the domains 
stay fixed in time. To this end we define the fixed Lagrangian domains = K 2 x (— m, 0) 
and fl + = IR 2 x (0, £) and assume that there exist invertible mappings 

7]° ± : n± ->■ fi±(0) (2.1) 

so that S = r/j_({x3 = 0}), 77+ ({X3 = £}) = {x% = £}, and r]°_ ({x 3 = — m}) = {x 3 = — m}. 
The first condition means that So is parameterized by the mapping 77° restricted to M 2 x {0}, 
and the latter two conditions mean that 77° map the fixed upper and lower boundaries into 
themselves. Define the flow maps r)± as the solution to 

dtV± (x,t) = u± {r]± (x, t),t) 
r](x, 0) = r]±(x). 

We think of the Eulerian coordinates as {y,t) with y = r](x,t), whereas we think of La- 
grangian coordinates as the fixed (x, t) E Q x R + ; this implies that ft±(t) = r]±(fl±,t) and 
that S(t) = 77 + ({x 3 = 0}, t), i.e. that the Eulerian domains of upper and lower fluids are the 
image of Q± under the mappings r]± and that the free interface is parameterized by r] + (-,t) 
restricted to M 2 x {0}. In order to switch back and forth from Lagrangian to Eulerian coor- 
dinates we assume that ry±(-, t) is invertible. Since the upper and lower fluids may slip across 
one another, we must introduce the slip map S± : IR 2 x M + — > R 2 x {0} C M 2 x (—m,£) 
defined by 

S-(x 1 ,x 2 ,t) = 7]Z 1 ( r n+(x 1 ,X2,0,t),t) (2.3) 

and S + (-,t) = S'Z 1 (-,t). The slip map £_ gives the particle in the lower fluid that is in 
contact with the particle of the upper fluid at x — (xi,X2, 0) on the contact surface at time 
t. 
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We define the Lagrangian unknowns 



v±(x,t)=u±(r}±(x,t),t) ^ 24 ^ 
q±(x,t) = p±(r)±(x,t),t), 

both of which are defined for (x,t) G Q± x IR + . Since the domains Q± are now fixed, we 
henceforth consolidate notation by writing rj, v, q to refer to r]±,v±, q± except when necessary 
to distinguish the two; when we write an equation for r),v,q we assume that the equation 
holds with the subscripts added on the domains £l±. Define the matrix A via A T = (Dr])^ 1 , 
where D is the derivative in x coordinates and superscript T denotes matrix transposition. 
Then in Lagrangian coordinates, the equations for v,q,r) are, writing dj = d/dxj, V = 
(<9i, d 2 , <9 3 ), and tr(-) for the matrix trace, 

8 t 7] = v 

d t q + qtr(ADv) = (2.5) 
d t v + AW(h(q)) = -gAVris- 

Since the boundary jump conditions in Eulerian coordinates are phrased in terms of 
jumps across the surface, the slip map must be employed in Lagrangian coordinates. The 
jump conditions in Lagrangian coordinates are 

(v + (x 1 ,x 2 ,0,t) -v-(S-(xi,x 2 ,t),t)) ■n(x 1 ,x 2 ,0,t) = 
P + (q + (x 1 ,x 2 ,0,t)) = P4Q4S-(.Xi,X2,t),t)) 

where we have written n — u or], i.e. 

djjj+ x d 2 r]+ 
\OiT] + X d 2 r] + \ 

for the normal to the surface £(i) = r] + ({xs = 0},t). Note that we could just as well have 
phrased the jump conditions in terms of the slip map S + and defined the surface and its 
normal vector in terms of r^_. Finally, we require 

v-(xi : x 2: —m,t) ■ e 3 = v+(xi,x 2 ,£,t) ■ e 3 = 0. (2.8) 

Note that since d t r) = v, 

e 3 • rj + (x 1: x 2 ,£,t) = e 3 • r]1(x 1 ,x 2: £) + / e 3 • v+(xi, x 2 , £, s)ds = £, (2.9) 

Jo 

which implies that r) + (xi,x 2 ,£,t) G {x 3 = £} for all t > 0, i.e. that the part of the upper 
fluid in contact with the fixed boundary {x 3 = £} never flows down from the boundary. It 
may, however, slip along the fixed boundary since we do not require v + (xi,x 2 ,£,t) ■ t{ = 
for i — 1,2. A similar result holds for at the lower fixed boundary {x 3 = — m}. 
In the subsequent analysis it will be convenient to employ the notation 

[/J=/ + I{Z3=0}-/-|{X3=0} (2-10) 

for the jump of a quantity / across the set {x 3 = 0}. 
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2.2 Linearization in Lagrangian coordinates 

In the steady-state case, the flow map is the identity mapping, 77 = Id, so that v = u 
and q = p. This means that the steady-state solution, p$, constructed above in Eulerian 
coordinates is also a steady state in Lagrangian coordinates. Now we want to linearize the 
equations around the steady-state solution v — 0, rj — Id, q = p , for which S_ = Id{ X3=0 } 
and A = I, the 3x3 identity matrix. The resulting linearized equations are 

{d t r] = v 
d t q + po div v = (2.11) 
p d t v + V(P'(po)<?) = -gqe 3 ~ 5 , PoV(e 3 • 77). 

Here we have rewritten the third equation by utilizing (11. 8p to write 

V(h'( Po )q) = -V(P'( Po )q) - ^^e 3 = -V(P'(p )q) + ^ (2.12) 
Po Po Po Po 

and then multiplying by po, which does not vanish on (—m,£). 
The jump conditions linearize to 

{v • e 3 ] = and [P'(p )?l = 0, (2.13) 

while the boundary conditions linearize to 

u_(ari, x 2 , -m, t) ■ e 3 = v + (x 1 , x 2 , £, t) ■ e 3 = 0. (2.14) 



3 Construction of a growing solution to (12.111 ) 



3.1 Growing mode ansatz 

We wish to construct a solution to the linearized equations ( 12. lip that has a growing H k 
norm for any k. We will construct such solutions via Fourier synthesis by first constructing 
a growing mode for a fixed spatial frequency. 

To begin, we assume a growing mode ansatz, i.e. let us assume that 

v(x, t) = w(x)e xt , q(x, t) = q(x)e xt , rj(x, t) = fj(x)e xt . (3.1) 

Here we assume that A G C\{0} is the same above and below the interface. A solution 
with 3ff(A) > corresponds to a growing mode. Plugging the ansatz into (12.111) . we get the 
equations 

Xfj = w 

\q + p divw = (3.2) 

Xp w + V(P'(Po)<?) = ~gqe 3 - gp V(e 3 ■ fj). 

Eliminating the unknowns fj and q by using the first and second equations, we arrive at the 
time-invariant system 

X 2 p w - V (P'(po)po div w) = gp Q div we 3 - gp Vw 3 . (3.3) 
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We may also eliminate q in the jump conditions to get 

M = and [P'(p )po div wj = 0. (3.4) 

Here we have used the fact that A ^ to remove it from the second jump condition. The 
boundary conditions are 

w 3 (x 1 ,x 2 , -m) = w 3 (x 1 ,x 2 , £) = 0. (3.5) 
3.2 Horizontal Fourier transformation 

Since the coefficients of the linear problem (13. 3p only depend on the x% variable, we are free 
to take the horizontal Fourier transform (11.201) . which we denote with either : or J 7 , to reduce 
to a system of ODEs in 23 for each fixed spatial frequency. 

We take the horizontal Fourier transform of wi, w 2 , W3 in (I3.3P and fix a spatial frequency 
£ = (61,6) G Define the new unknowns <p(x 3 ) = £2, £3), #(^3) = iwjfc, £2, £3), 

and ^(#3) = u>3(£i, £2, ^3) so that 

JF(divw) = fi¥> + &0 + (3-6) 
where ' = d/dx 3 . Then for <p, 8, ip and A = A(£) we arrive at the following system of ODEs. 

(A 2 po + £ 2 P'(po)po) <P = -6 (P'(Po W - 9Po4>) - £i6P'(pq)pq0 (3.7) 

(a 2 P o + e 2 p'(Po)po) = -6 (p'(po w - ^po^) - eie 2 p'(po)po<p (3.8) 

- (P'( P o W)' + *W = 6 ((P'(poW)' + <?Po<p) + 6 {(P'(po)p 8)' + SPo#) (3.9) 

We can reduce the complexity of the problem by removing the component 9. To do 
this, note that if p,0,ip solve the above equations for £1,^2 and A, then for any rotation 
operator R 6 5*0(2), (<p, 6) := R(<p,$) solve the same equations for (£1,^2) : = P(£i?£2) with 
ip, A unchanged. So, by choosing an appropriate rotation, we may assume without loss of 
generality that £2 = and £1 = |£| > 0. In this setting 9 solves 

\ 2 p 9 = 0, (3.10) 

and since po > on (— m, £) we have that 6 = 0. This reduces to the pair of equations 

-aw = iei (p'(po)po^ + lei v)) - 9 lei p^ , 3 n) 

-A 2 p ^ = -(p / (po)po(^ + lei <p))' - 9 lei pov? 

with the jump conditions 

M = and [P'(po)po(|£| + = (3.12) 

and the boundary conditions ip(—m) = ip(£) =0. 

Since A G C, <p, t/> can be complex valued, but we may reduce to the case of real valued 
functions as follows. Multiply the first equation in (13. HI) by <p and the second by ip, add the 
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resulting equations, and integrate over (— m, £). An integration by parts and an application 
of the boundary and jump conditions shows that 



- tt / MM* + M 2 ) = \ f P'(Po)Po W + ICI ¥>f - 2g ICI Pb»W). (3-13) 

From this we see that —A 2 G K, which implies that 3?(y), 3ft(V ; ) are also solutions, so we may 
restrict ourselves to finding real valued solutions. Moreover, we know that either A G R if 
—A 2 < or A G iR if —A 2 > 0. Since we want A > for a growing mode solution, we are 
interested only in finding solutions with —A 2 := p < 0. For a given |£|, the smallest value of 
—A 2 = p (and hence the largest growing mode if p < 0) can then be found by a variational 
principle: 

M = A*(fcl) 

1 



1 * 



= mf | - y p / (po)po(^ / + lei - 2^ lei p^> 

We will use this in the next section to find solutions. 



2 



Po (^ + ^) = l . (3.14) 



3.3 Constrained minimization formulation 

In this section we will produce a solution to ( 13. lip with |£| > by utilizing variational 
methods. In order to understand A in a variational framework we consider the two energies 

1 r e 

E{ip, il>) = - / P'(p )paW + 1^1 ^) 2 - 2(7Po If | </><P (3.15) 

^ J-m 

and 

JM) = \( Po(^ 2 + ^ 2 ), (3.16) 

which are both well-defined on the space L 2 ((—m,£)) x Hq((— m, £)). In this section we 
will write Hl((—m,£) for the usual Sobolev space of functions on (—m,£) that vanish at 
the endpoints, i.e. for our variational methods we will not use the piecewise Sobolev space 
defined in the introduction. Consider the set 

A={M) G L\(-m,£)) x Hl{{-m,£)) | J(^) = 1}. (3.17) 

We will generate solutions by minimizing E over the collection A. In order for these 
minimizers to give rise to a solution to (13. lip , it must hold that 

-A 2 = inf E((p,ip)<0. (3.18) 

Throughout this section we assume that |£| > 0, and we will construct a non-trivial solution 
with A = A(|£|) > 0. To this end, we begin by showing that the infimum is negative for each 
|£|>0. 
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Lemma 3.1. It holds that m£{E(ip,ip) \ (ip,ip) G A} < 0. 

Proof. Since both E and J are homogeneous of degree 2 it suffices to show that 

E(cp,ip) 



inf 



< 0, 



(3.19) 



but since J is positive definite, we may reduce to constructing any pair (<p, ip) G L 2 x ifg 
such that E(ip,ip) < 0. We will further assume that <p = —ip' / |£| so that the first term in 
E((p, ip) vanishes. We must then construct ip G Hq so that 



E(iP):=E(-iP'/\Z\,iP) = / gPoW>'<0- (3.20) 
We employ the identity ipip' = {ip 2 ) 1 /2 and an integration by parts to write 



£Po^ 2 



2 ■/!"> 



^Po^ 2 + 



£Po^ 2 



giP 2 (Q) 



M + 



Po 



^'(Po) 



ip 2 . (3.21) 



Notice that by f)1.14p . [p ] = Po — Po > so that the right hand side is not positive definite. 
For a > 2 we define the function ip a G Hq{(— m,£)) according to 



1pa(x 3 ) 



V+p;(po(x 3 )) 

' PqPL( P q(x 3 )) 
^PL(Po)po(x 3 ) 



) 1/2 (l-f)° /2 , x 3 G(0,£) 



(3.22) 



According to (11. lip . ip a is continuous across x% = and -?/> a (0) = 1. A simple calculation 
then shows that 

1 / mpo 



Po /2 



P'(p 



+ 



1 + a \PL( Po ) ' P|(p^ 



which implies that 



[Po] + 



mp 



+ 



(3.23) 



(3.24) 



a + HPi(Po) P'M 
Then for a sufficiently large we have that E{ip a ) < 0, thereby proving the result. □ 
The key to the argument presented in Lemma [3. II was constructing a pair (<p,ip) so that 



-9 IPo] ,2 



r(o) < o, 



(3.25) 



which in particular required that ^(0) ^ 0. We can show that this property is satisfied by 
any (<p, ip) G A so that E((p, ip) < 0. 
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Lemma 3.2. Suppose that (<p,ip) £ A. satisfy E(ip,ip) < 0. Then ip(0) 7^ 0. 
Proof. A completion of the square allows us to write 

P'(p )PoW + \Z\<p) 2 -2gpo\Z\^ 

y/ p (Po) J P (Po) 

Integrating by parts as in (13.211) . we know that 

2 9t>oW - Sr^-y = -9 IM ^ 2 (o). (3.27) 

Combining these equalities, we can rewrite the energy as 

e{<p, vo = -f m ^ 2 (o) + ~ f_ m (vn^ioW + if i & - -jj^A ■ ( 3 - 28 ) 

From the non-negativity of the integrals, we deduce that if E(ip,ip) < 0, then ip(0) 7^ 0. □ 

With Lemma 13.11 in hand, we way apply the direct methods to deduce the existence of a 
minimizer of E in A. 

Proposition 3.3. E achieves its infimum on A. 

Proof. First note that by employing the identity — lab = (a — b) 2 — (a 2 +b 2 ) and the constraint 
J(ip, ip) — 1 we may rewrite 

E(cp, ip) = -g |f | + ~ / P'( Po )poW + If | V) 2 + 9 If I Po(v ~ > -9 If I • (3.29) 

^ J— m 

This shows that E is bounded below on A. Let (<p n , ip n ) G A be a minimizing sequence. Then 
<y? n is bounded in L 2 ((— m, £)) and ip n is bounded in i?o((— m, £)), so up to the extraction of 
a subsequence ip n — </? weakly in L 2 , ■?/>„ — ^ ^ weakly in ifg, and ^ n — >■ ip strongly in L 2 . Via 
weak lower semi-continuity and the strong L 2 convergence ip n — > ip, we then have that 

E(cp, ip) < lim inf E(<p n , ip n ) = inf (3.30) 

All that remains is to show that (ip, ip) G *4. 

Again by lower semicontinuity, we know that J(<p,ip) < 1. Suppose by way of contradic- 
tion that J((p,ip) < 1. By the homogeneity of J we may find a > 1 so that J(aip,aip) = 1, 
i.e. we may scale up (ip, ip) so that (ay?, cap) G A By Lemma I3TT1 we know that inf^E < 0, 
and from this we deduce that 

E(a<p, oiip) = a 2 E(y, ip) < a 2 inf E < inf E, (3.31) 
which is a contradiction since (a<p, aip) G A. Hence J(ip, ip) = 1 so that (<p, ip) G A. □ 
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Remark 3.4. By inequality f!3.29j) we know that —A 2 = E{ip,ip) > — g |£| and hence A < 



g\t\- 

The next result provides a lower bound for A for large |£|, showing that A(|£|) — > oo as 

Lemma 3.5. There exist constants Cq, Ci, Ci > depending on p^ , P±, g,m, £ such that the 
eigenvalue A = A(|£|) satisfies 

A 2 > d |d - C 2 /or |d > C . (3.32) 
Proof. For a > 2 we consider the function ^ a defined by 

(l-f) a/2 , x 3 G[0,£) 
(l + f) a/2 , x 3 e(-m,0). 



Then according to (I3T2H and (EEZfll 

lei > < "A + t4^ ( 3 - 34 ) 

1 + a 

for two constants Al, A 2 > depending on p , P±, o, m, £ but not on a or |£|. On the other 
hand, since po is bounded below and above, straightforward calculations yield the bounds 

J(-<l \i\M = \ £ po (tf + <g£) < ^ + Jj, + * + (3 ' 35) 

and 

lei > ( 3 - 36 ) 

1 + a 

for constants A 3 , A4, A 5 , Aq, A 7 > depending on the same parameters but not a or |£|. 
From this we see that for a = |£| > Co sufficiently large, there are positive constants Ci, C 2 
so that 

|=#f|[44<-C lk | + C , (3.37) 

Since —A 2 = inf E/J, the result follows. □ 

We now show that the minimizers of Proposition 13.31 satisfy Euler-Langrange equations 
equivalent to (13. lip . 

Proposition 3.6. Let (<£>, VO G A be the minimizers of E constructed in Proposition \ 3.3[ 
Let p := E{ip,ip) < 0. Then (<p,ijj) satisfy 

'pp ip = iei (p'(po)po(^ + lei <p)) - 9 lei p v> f3 38) 

= -(p'(po)po(^ + iei v))' - 9 iei po^ 

a/ong i/ie jitmp conditions 

M = and [P'(p )po(|el V + = (3.39) 

and £/ie boundary conditions ip(—m) = ip(£) = 0. Moreover, the solutions are smooth when 
restricted to either (— m, 0) or (0, £). 



14 



Proof. Fix (ifo^o) e L 2 ((-m,£)) x H%((-m,£)). Define 

j(t, s) = J{<p + tip + s(f,ip + tip + sip) (3.40) 
and note that j(0, 0) = 1. Moreover j is smooth and 

§(0,0) = / p [w + VW>] (3.41) 



ft 
and 

§(0,0) = / po(^ + ^ 2 ) = 2. (3.42) 
ds J_ m 

So, by the inverse function theorem, we can solve for s = o~(t) in a neighborhood of as a 
C l function so that cr(0) = and j(t,a(t)) = 1. We may differentiate the last equation to 
find 

^(0,0) + j|(0,0y(0) = 0, (3.43) 

and hence that 

Idj' 1 " 



o'(O) = ~2^(°' °) = "a / Po ^ + ^ • (3 - 44) 



Since (if, if>) are minimizers over A, we then have 

£(y? + tfo + a(t)ip, ip + tip + tr{t)ij}), (3.45) 



d 

= -r 

dt 



t=0 



which implies that 



o = / p'( Po )poW + iei v)(^o + ^' w + lei + lei ^(om 

J —m 

- 9 If I Po(^(y?o + ct'(0)^) + V (4j + tT'(O)V)). (3.46) 
Rearranging and plugging in the value of cr'(0), we may rewrite this equation as 

P'(Po)PqW + Ifl ¥>)(V£ + If I <Po) - 2 If I Po(^o + <?A) =H f Po [PW + <M (3.47) 

J —m 

where the Lagrange multiplier (eigenvalue) is \i = E((p,ip). 

By making variations with <p , ip compactly supported in either (— m, 0) or (0, £), we find 
that ip and ip satisfy the equations (I3.38P in a weak sense in (— m, 0) and (0, £). Standard 
bootstrapping arguments then show that (y?,V0 are i n H k ((—m,0)) (resp. H k ((0,£))) for all 
k > when restricted to (— m, 0) (resp. (0,£)), and hence the functions are smooth when 
restricted to either interval. This implies that the equations are also classically satisfied 
on (— m, 0) and (0,£). Since ((p,if>) € H 2 , the traces of the functions and their derivatives 
are well-defined at the endpoints x% = —m,0,£, and it remains to show that the jump 
conditions are satisfied at x% = and the boundary conditions satisfied at x^ = —m, £. 
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Making variations with respect to arbitrary ip , ip Q £ C%°((—m,£)), we find that the jump 
condition 

[P'(p )po(|£k + VOJ = 0. (3.48) 

must be satisfied. Note that the conditions [p0] = and ip(—m) = tp(£) = are satisfied 
trivially since ip G H%((-m, £)) ^ C °' 1/2 ((-m, £)). □ 

We may now deduce the existence of a growing mode solution. 

Theorem 3.7. For any |£| > there exists a solution ((f,ip) with A = A(|£|) > to (13. lip 

satisfying the appropriate jump and boundary conditions. For these solutions ip(0) ^ 0, and 
the solutions are smooth when restricted to (— m, 0) or (0,1). 

Proof. Let A 2 (£) = —y, = —E(<p, ip) > 0, where (<p, tjj) are the minimizers of E over the set A 
constructed in Proposition 13.31 Then Proposition 13.61 shows that (</?, iji) solve the equations 
( 13~TT|) with A = y/^JI > 0. Lemma EZ3 implies that ^(0) ^ 0. □ 

The next result provides an estimate for the H k norm of the solutions ((p,ip) with |£| 
varying, which will be useful in the next section when such solutions are integrated in a 
Fourier synthesis. To emphasize the dependence on |£| we will write these solutions as 

<p(\Z\) = <p(\Z\,x 3 ) andV(i£|)=V(l£i,Z3). 

Lemma 3.8. Let (p(\£\), V'(lCl) ^ e ^ e solutions to (13. 1 1 j) constructed in Theorem \3. 7[ Lei 
Co, Ci, C*2 > 6e £/ie constants from Lemma \3. 51 Fzx < -Ri < oo so i/ia£ 

i?i > max{C , 2C 2 /C X }. (3.49) 

JTien /or |£| > R lt /or eac/i k > t/iere exists a constant > depending on P±, g, m, £ 
so t/iat 

!Mi£i)iu (( - OTl o)) + ii^(i^i)L fc ((- m ,o)) + i^(iei)iifl*( ( o^ + wmDWHHm) < ^Eitf- 

i=o 

(3.50) 

v4feo, i/iere exists a B > depending on the same parameters so that for any |£| > 

vV(l£l)+^ 2 (l£l) r2 „ m > Bb- (3-51) 

L^ix— m,l)) 

Proof. We begin with the proof of (I3.50p . For simplicity we will prove an estimate of the 
H k norms only on the interval (0,£). A bound on (— m, 0) follows similarly, and the result 
follows by adding the two. First note that the choice of Ri, when combined with Lemma 
13.51 and the bound (I3.29P implies that 

-0|fl<A*<tf a -tfi|e|<-y|£|. (3-52) 

Recall also that p is smooth on each interval (0,£) and (— m, 0) and bounded above and 
below. Throughout the proof we will let C > denote a generic constant depending on the 
appropriate parameters. 
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We proceed by induction on k. For k = the fact that (<£>(|£|), ^(ICD) £ A implies that 
there is a constant A > depending on the various parameters so that 

ym)\\ L H m ) + mm L2m) )< A °- (^) 

Suppose now that the bound holds some k > 0, i.e. 

k 

+ IIV>(l£i)IU( ( o,*)) <^E Itf • ( 3 - 54 ) 

Define |) := P;(po)Po(^'(l£l) + If I ¥>(ICI)), where ' = 9 X3 . Then ([331 implies that 

'^{\i\) = ^\iV 1 pM\i\) + 9PMi\) (3 55) 

<J{\t\) = -mm\)-g\t\fw{\Z\)- 

These equations and (13.521) then imply that 

H k ((o,e)) <C(M\m H H(0A) + \\^\\HHm)) ( 3 - 56 ) 

and 

ll^'Cie|)II^ C<: o^)) < Cl^l (!I^Cie|)Il^ fcC<: o^)) l|-0(|e|)||^ feao ^))) (3-57) 
so that ||w(|f|) ||jjfc+i((o,«) — C^jio IfP- But then the definition of w(|£|) implies that 



HH(o,e)) 



so that H^(l^l) Hh'=+ 1 ((o £)) — C^^io If P- Returning to the first equation in (13.551) . we see 
that 



i^(iei)L, +1 ( ( o ) ,))<^ii^(iei)iu +1 ((o ) .)) + f 



(3.59) 

fe+i 



Invoking the bound (13.521) again, we deduce that ||v ? (I£I)IIh'=+i((o£)) — CE 3 =o IfP '• Hence 

fe+i 

<A fc+1 ^|eP (3.60) 
j=o 

for some constant A^+i > depending on the parameters, i.e. the bound holds for k + 1. 
By induction, the bound holds for all k > 0. 

To prove (13.511) we again utilize the fact that for any |£| > we have (<p(|f|), VKIfD) £ -4- 
Since p is bounded above and below, the bound follows. 

□ 

A solution to (13. lip gives rise to a solution to the system of equations (13 . 7[) ( T379T) for the 
growing mode velocity, w, as well. 
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Corollary 3.9. For any £ = (£1, £ 2 ) 7^ (0, 0) there exists a solution tp = x 3 ) ; 6 1 = x 3 ), 
V* = V'dCl 5^3) anc ^ Ml£l) > to (j3.7p - (l3.9|) satisfying the appropriate jump and boundary 
conditions so that i[)(0) 7^ 0; £/ie solutions are smooth when restricted to (— m, 0) or (0, £). 
TTje solutions are equivariant in £ in £/ie sense t/iat if R & SO (2) zs a rotation operator, then 

<p(R£,x 3 )\ [R u R12 0\ [<pfox 3 )\ 

e(RC,x 3 ) \ = \R 21 R22 flfca*) . (3.61) 

^m,x 3 )j \ 01/ \v(e^3)/ 

Proof. We may find a rotation operator i? G 50(2) so that _R£ = (|£| ,0). Define 

(¥>(£, x 3 ), 0(£, x 3 )) = ^Miei , *s), 0) (3.62) 

and i[)(£,x 3 ) = ) x 3)> where y?(|£|) and ^(ICI) are the solutions from Theorem 13.71 This 
gives a solution to (j3.7j) — (j3.9j) . The equivariance in £ follows from the definition. □ 



3.4 Fourier synthesis 

In this section we will use Fourier synthesis to build growing solutions to (12.111) out of 
the solutions constructed in the previous section (Corollary 13 .9[) for fixed spatial frequency 
(Gl 2 . The solutions will be constructed to grow in the piecewise Sobolev space of order k, 
H k } defined by ffL22|) . 

Theorem 3.10. Let Ri < R 2 < R 3 < 00, where Ri > is the constant from Lemma 
\3.8l Let f G 0^°(M 2 ) be a real-valued function so that /(£) = /(|£|) and supp(/) C 
B{0,R 3 )\B(0,R 2 ). For^eR 2 define 

^3) = x 3 )ei - i8(£, x 3 )e 2 + ^(£> x 3 )e 3 , (3.63) 

where (p,9,ip are the solutions provided by Corollary \3.9i . Writing x' ■ £ — Xi£i + x 2 C, 2 , we 
define 

r,(x,t) = -± 1 f m^,^y mt e ix '^, (3.64) 

47T ./ TO 2 

47T 



v(x,t) = -^- I \mf(tWt,^ m <P*dt, (3-65) 



and 

Po(x 3 ) 



q(x,t) 



47T 5 



/ f(O^M^x 3 )+^,x 3 ) + d x ^,x 3 ))e x ^ t e^'<d^ (3.66) 

Then rj, v, q are real-valued solutions to the linearized equations (12. lip along with the corre- 
sponding jump and boundary conditions. The solutions are equivariant in the sense that if 
R G 50(3) is a rotation that keeps the vector e 3 fixed, then 

rj(Rx, t) = Rr](x, t),v(Rx, t) = Rv(x, t), and q(Rx, t) = q(x, t). (3.67) 

For every k G N we have the estimate 

\\v(0)\\ Hk + M0)\\ B " + ll<K0)ll*» < (jT a (l + le| 2 ) fc+1 1/(0! 2 ^) 1/2 < 00 (3.68) 



for a constant Ck > depending on the parameters p^, P±,g, m, I; moreover, for every t > 
we have rj(t),v(t),q(t) e H k (Q) and 

\\ V (0)\\h- < Wv(t)\\ H * < e tVm \\v(0)\\ Hk 
sVchR^c; \\ v (Q)\\ Hk < \\v(t)\\ Hk < \\v(0)\\ Hk (3.69) 



e 



WC\R 2 -C 2 



ll9(0)||^<||9(t)||^<e*^||?(0)|U 



where C\, C 2 > are the constants from Lemma l375\ and \JC1R2 — C2 > 0. 
Proof. For each fixed (el 2 , 

v(x,t) = f(£)*(t,x 3 )eWMe«< (3.70) 

v(x, t) = \(\Z\)m)w(t, x 3 )e x ^ t e i ^, (3.71) 

and 

q(x,t) = -po{x 3 )f(€)($MZ,x 3 ) +60(£,x 3 ) + d 3 ^,x 3 ))e x( -^ t e ix ^ (3.72) 
gives a solution to (12. lip . Since supp(/) C B(0, R 3 )\B(0, R 2 ), Lemma [3781 implies that 

sup ||fl£ «)(£,. )|L < 00 for all fceN. (3.73) 



5esupp(/) 



Also, A (|(|) < \J g These bounds imply that the Fourier synthesis of these solutions given 
by ( 13.64p ~ 113.66p is also a solution to (12.1ip . The Fourier synthesis is real- valued because / 
is real- valued and radial and because of the equivariance in £ given in Corollary 13.91 This 
equivariance in £ also implies the equivariance of rj,v,q written in 1 13. 6 7]) . 

The bound ( I3.68P follows by applying Lemma 13.81 with arbitrary k > and utilizing the 
fact that / is compactly supported. According to ( 13. 52ft 



^fgRz > vW\ > A(|£|) > s/Gy |f | - C 2 > s/CJii - C 2 > 0, (3.74) 
which then yields the bounds (13.691) . □ 
Remark 3.11. It holds that 

^(^,^,0,0) = ---/ f(t)4>(£,0)e"*d£ (3.75) 

is the vertical component of the initial linearized flow map at the interface between the two 
fluids. Since ?/>(£, 0) 7^ for any choice of a nonzero f in general gives rise to a nonzero 
r] 3 (x 1 ,x 2 ,0,0). 



4 Analysis of the linear problem 

4.1 Estimates for band-limited solutions 

Suppose that rj,v,q are real-valued solutions to ( 12. lip along with the corresponding jump 
and boundary conditions (of course, by linearity, we may also handle complex solutions by 
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taking the real and complex parts and proceeding with an analysis of each part). Further 
suppose that the solutions are band-limited at radius R > 0, i.e. that 

|J supp(|r)(-,x 3 )| + \v(;x 3 ))\ + \q(-,x 3 )\) C B(0,R) (4.1) 

X3£(— m,t) 

where v denotes the horizontal Fourier transform defined by (ll.20p . The importance of the 
band-limited assumption lies in the fact that the growing mode exponent remains bounded 
for |£| < R. We henceforth denote 

~ (4.2) 



A(R) : = sup A(|£|) < y/gR < oo. 

0<\£,\<R 

We will derive estimates for band-limited solutions in terms of A(R). 

It will be convenient to work with a second-order formulation of the equations. To arrive 
at this, we differentiate the third equation in (12. lip in time and eliminate the q and rj terms 
using the first and second equations. This yields the equation 

poduv - V(P'(p )po div v) + gp Q Vv 3 - gp div ve 3 = 0. (4.3) 

coupled to the jump conditions 

ld t v 3 j = and [P'(p )po div vj = (4.4) 

and the boundary condition d t v 3 (xi,X2, —m,t) = d t v 3 (xi,X2,i,t) = 0. The band limited 
assumption implies that supp(t>(-, x 3 )) C B(0,R) for all x 3 G (—m,£). The initial data for 
d t v (0) is given in terms of the initial data g(0) and r/(0) via the third linear equation, i.e. 
d t v (0) satisfies 

p d t v(0) = -V(P'(p )?(0)) - gq(0)e 3 - ^p Vr ?3 (0). (4.5) 
Our first result gives an evolution equation for an energy associated to such solutions. 
Lemma 4.1. For solutions to (14. 3|) it holds that 



p0 -\d t v\ 2 + P ' ip0)p0 



div v — 



P'(pc 



-v 3 



gM I . ,2 



(4.6) 



Proof. Recall that f2 + = M 2 x (0, £). Take the dot product of (14. 3 p with d t v(t) and integrate 
over Q + . After integrating by parts and utilizing (II .8p . we get 

/ p Q d t v ■ d tt v + P'(Po)po(div u)(div d t v) - gp (v 3 div d t v + d t v 3 div v) + ^., P °. v 3 d t v 3 
Ja+ P'(po) 



gpoV 3 d t v 3 



P;(p + )p + div^ 3 . (4.7) 



We may then pull the time derivatives outside the first and second integrals to arrive at the 
equality 



Po- 



\d t v\ 2 , P'(Po)po 



divf 



g 



P'(Po) 



v 3 



gpo 



\v 3 \ 



P' + {pt)pt div vd t v 3 . (4.8) 
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A similar result holds on Q_ = M? x (— m, 0) with the opposite sign on the right hand side. 
Adding the two together yields 



d t / Po 



\d t v\ 2 , P'(Po)Po 



div v — 



P'(Po) 



v 3 



dt [ 9 [poj 
Jr 2 



\ v 3\ 



lP'(p )p divvd t v 3 }. (4.9) 



The jump conditions imply that the second term on the right hand side of this equation 
vanishes, and the result follows. □ 

The next result allows us to estimate the energy in terms of A(R). 

Lemma 4.2. Let v G H 1 ^) be band-limited at radius R > and satisfy the boundary 
conditions v 3 (xi, x 2 , —m, t) = v 3 (xi,x 2 , £, t) — 0. Then 



g{po 



v 3 



P'{po)Po 



divt> 



P'iPo 



-v 3 



< 



A 2 (R) 



Po \v\ 



(4.10) 



where A(R) is defined by (14.21) . 

Proof. Take the horizontal Fourier transform ( ll.20p and apply (ll.2ip to see that 



An 2 



glpo 



N -4vr 2 



9 [Po] ,2 
- \v 3 \ di 



P'(po)po 



divv 



P'iPc 



-v 3 



P'(po)po 



i£ivi + i&h + 9 X3 v 3 - v 3 

P(po) 



9lPo\ 



V 3 



P\po)Po 



i^ivi + i&v 2 + d X3 v 3 



P'iPo 



-v 3 



d£dx 3 
dx 3 ) d£. (4.11) 



Consider now the last integrand for fixed (£1,^2) = £ 7^ 0, writing (p(x 3 ) = u>i(£i, £2, x 3 ) 
9(x 3 ) = iv 2 (£i,&,x 3 ), i>(x 3 ) = v 3 (£, h £ 2 ,x 3 ). That is, define 



where ' = d X3 . By splitting 



P'iPo)Po 



P'iPo) 



dxi 



Z(<p, 9, if>; = Z(&<p, W, 3fty>; f) + Z(Q(p, 30, 3^>; f) 



(4.12) 



(4.13) 



we may reduce to bounding Z when (p, 9, ip are real-valued functions, and then apply the 
bound to the real and imaginary parts of (p, 9, ip. 

The expression for Z is invariant under simultaneous rotations of (£1,^2) and {(p,0), so 
without loss of generality we may assume that (£1, £ 2 ) = (|£| , 0) with |£| > 0. Then according 
to (jS^HD , Z(<p,6,ip;€) = -Efaip) and hence 



Zi<p,9,i>;0< 



A 2 (ie|) 



Po(<fi 2 + ^ 



(4.14) 
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where A(|£|) is the growing-mode exponent constructed in the previous section. Translating 
the inequality back to the original notation for fixed £, we find 



— ^— N - 



i£iVi + i&h + d 3 v 3 - — 



P'(Po) 



v 3 



dx 3 < 



A 2 (|CI) 



Po \v\ 



(4-15) 

Integrating each side of this inequality over £ 6 IR 2 , employing the fact that A 2 (|£|) < A 2 (R) 
on the support of v, and using (11.211) then proves the result. □ 

We may now derive growth estimates in terms of the initial data and the value of A(R). 

Proposition 4.3. Let v be a solution to (14. 3 p along with the corresponding jump and bound- 
ary conditions that is also band-limited at radius R > 0. Then 



lia (n) + \\dtv(t)\\ L2(n) 

< Ce 2A ^ (\\v(0)\\Un) + \\dtv(0)\\l m + \\^ v(0)\\ 2 L2(n) ) (4.16) 
for a constant < C = C(p^, P±, g, A(R), m, t). 

Proof. Integrate the result of Lemma 14.11 in time from to t to find that 



Po ' „' <A + 



where 



.4 



Po ^ 1- 



^s(t)| 2 - 

^'(Po)po 



^'(Po)Po 



divv(t) 



P'iPc 



divv(0) 



P'(Po) 



«3(0) 



We may then apply Lemma 14.21 to get the inequality 

H0 2 2 

|2 r ixi2 



po K*)r 



(4.17) 



(4.18) 



(4.19) 



Defining the weighted L 2 norm by 
inequality as 

\\\dXt)\\ 2 <A + ^\Wt)\\ 2 - 

An application of Cauchy's inequality shows that 

A(R)d t \\v(t)\\ 2 = A(R)2(d t v(t),v(t)) < A 2 (R) \\v(t)f + IIM*)I| 2 



Jo Po \ f\ ? we ma Y compactly rewrite the previous 



(4.20) 



(4.21) 



where (•, •) is the weighted L 2 inner-product associated with the weighted L 2 norm. We may 
combine this inequality with (14.201) to derive the differential inequality 



d t \\v 



2 < B + 2A(R) \\v(t)\\ 2 



for B = 2A/A(R). Gronwall's inequality then shows that 



Kt)|| 2 <e 2A ^|M0)|| 2 + 



B 

2A(R) { 



2A{R)t _ 



1) 



(4.22) 



(4.23) 
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for all t > 0. To derive the corresponding bound for ||9 t u(t)|| 2 we return to (I4.20p and plug 
in (TPS]) to see that 

\ \\dMt)f < A + ^ (e 2A ^ |K0) || 2 + ^y(e 2A(i?)t - 1)) • (4.24) 

The result follows by noting that the weighted L 2 norm ||-|| is equivalent to the standard L 2 
norm ||-|| L 2(q) and that 

A < C (\\v(0)\\ 2 L2(n) + ||^(0)||* 2(n) + ||di™(0)|| 2 2(fi) ) (4.25) 
for a constant C > depending on the parameters pg, P±, g,m,£. □ 

4.2 Uniqueness 

In order to employ Proposition 14.31 in a proof of uniqueness for solutions to the linearized 
equations, we must first construct a horizontal spatial frequency projection operator. To 
this end let $ G C C °°(M 2 ) be so that < $ < 1, supp($) C 5(0,1), and $(x) = 1 for 
x G 5(0, 1/2). For R > let $ K be the function defined by ® R (x) = <&{x/R). We define the 
projection operator Pr via 

P R f = 7- 1 {<$> R Ff) (4.26) 

for / G L 2 (fl), where J 7 = A denotes the horizontal Fourier transform in Xi,x% defined by 
fll.20p . It is a simple matter to see that P R satisfies the following. 

1. P R f is band- limited at radius R. 

2. P R is a bounded linear operator on H k (Q) for all k > 0. 

3. P R commutes with partial differentiation and multiplication by functions depending 
only on £3. 

4. P R f = for all R > if and and only if / = 0. 

With this operator in hand we may prove the uniqueness result. 
Theorem 4.4. Solutions to ( 12. lip are unique. 

Proof. It suffices to show that solutions to (12.111) with initial data remain for t > 0. 
Suppose then that rj, v, q are solutions with vanishing initial data. Fix R > and define 
Vr = PrV, v R = Pr v i Qr = PrQ- The properties of P R show that r) R ,v R ,q R are also 
solutions to (12. lip but that they are band-limited at radius R. Turning to the second order 
formulation, we find that v R is a solution to (14.31) with initial data v R (0) = d t v R (0) = 0. We 
may then apply Proposition 14.31 to deduce that 

IM*) 11^(0) = \\dtV R (t)\\ L2{n) = for all t > 0, (4.27) 

which implies that r] R (t),v R (t), and q R (t) all vanish for t > 0. Since R was arbitrary it must 
hold that T](t), v(t) and q{t) also vanish for t > 0. □ 
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4.3 Discontinuous dependence on the initial data 

The solutions to the linear problem (12. lip constructed in Theorem 13.101 are sufficiently patho- 
logical to give rise to a result showing that the solutions depend discontinuously on the initial 
data. Then, in spite of the previous uniqueness result, we get that the linear problem is ill- 
posed in the sense of Hadamard. 

Theorem 4.5. The linear problem (12. lip with the corresponding jump and boundary con- 
ditions is ill-posed in the sense of Hadamard in H k (Q) for every k. More precisely, for any 
k,j G N with j > k and for any T > and a > there exists a sequence of solutions 
{( r ln,v n ,q n )}'^' =l to (12. lip , satisfying the corresponding jump and boundary conditions, so 
that 

IM0)|U + K(0)|U + |M0)|| ffJ < i (4.28) 

but 

\K{t)\\ H k > \\Vn{t)\\ H u > a for all t > T . (4.29) 

Proof. Fix j > k > 0, a > 0, T > and let Cj, R±, B , Ci, C 2 > be the constants 
from Theorem I3.10[ Lemma 13.81 and Lemma 13.51 respectively. For each n G N let R(n) be 
sufficiently large so that R(n) > R±, ^JC\R{n) — C 2 > 1, and 



exp(2T 0V /^H-C 2 ) q 

(1 + (R(n) + ~ an Bl ^ ] 

Choose /„ G C™(R 2 ) so that supp(/ n ) C B(0,R(n) + 1)\S(0, R{n)), f n is real-valued and 
radial, and 



We may now apply Theorem 13.101 with f n , R 2 = R(n), and R 3 = R{n) + 1 to find 
Vn,v n ,q n that solve (12.111) and the corresponding jump and boundary conditions so that 
Vn{t) , v n (t) , q n (t) G fP(f2) for all t > 0. By (13.681) and the choice of satisfying (14.311) . we 
have that (j4~28l) holds for all n. 

We may then estimate 



lk(T )|| 2 H , > / (l + |e| 2 ) fc |/n(0| 2 e 2toA(l?l) |l^,-)lli 2( _ m/) ^ 

> M% [ (1 + If | 2 )' +1 1/„(? )| 2 = a 2 . 

Here the first bound is trivial, the second follows since supp(/ n ) C B(0,R(n) + 1) and 
A(|£|) > a/ CiR(n) — C 2 , and the third follows from the choice of R(n) and the lower bound 
(I3.5ip . Since A(|£|) > a/ CiR(n) — C 2 > 1 on the support of f n , we also know that 

\\v n (t)\\ 2 Hh > \\Vn(t)\\ 2 H , > \\Vn(T )\\ 2 H k for * > T , (4.33) 
from which we deduce (I4.29p . □ 
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5 Ill-posedness for the non-linear problem 

Recall that the steady state solution to (12. 5p is given by v = 0, r\ = r]^ 1 = Id, q = p with 
A = I and >S_ = S + = Id{ X3=0 }. We will now rephrase the non-linear equations (12. 5 j) in a 
perturbation formulation around the steady state. Let 

rj = Id + fj, "if 1 = Id - C, v = + v, q = po + a, A = I - B, (5.1) 

where 

oo 

B T = ^(-ir 1 ( J Dr/)". (5.2) 

n=l 

In order to deal with the term h(q) = h(p + er) we introduce the Taylor expansion 

h(p + a) = h(p ) + h'(p )a + n (5.3) 
where the remainder term is defined by 

IZ{x,t) = / (a(x,t)-z)h"(p (x)+z)dz = / (p (x)+a(x,t)-z)h"(z)dz. (5.4) 



Then the evolution equations (I2.5P can be rewritten for fj, v, o as 

dtfj = v 

d t a+(p + cr)(dwv-tr(BDv))=0 (5.5) 
d t v + (I- B)V(h'( Po )a + gfj 3 +K)=0. 

We require the compatibility between ( and fj given by 

C = 77o(/d-C) (5.6) 

The jump conditions across the interface are 



(v + (xi,x 2 ,0,i) - v_(S^(xi,x 2 ,t),t)) ■ n(x 1 ,x 2 ,0,t) = 
P+(Po + a+(xi,x 2 ,0,t)) = P-(p + a-(S-(x 1 ,x 2 ,t),t)) 



(5.7) 



where the slip map (I2.3P is rewritten as 

S- = (Id R 2 - £_) o (Jd R 2 + r/+) = /d R 2 + - C_ o (Jd R2 + r] + ). (5.8) 

Finally, we require the boundary condition 

V-(xi,X2, —m,t) ■ e 3 = v + (xi,x 2 ,£,t) ■ e 3 = 0. (5.9) 

We collectively refer to the evolution, jump, and boundary equations (I5.5l) - (l5.9p as "the 
perturbed problem." To shorten notation, for k > we define 

||(t;,i;,<7)(*)|| h * = ||r7(t)||^ + \\v(t)\\ Hh + \\a(t)\\ Hk . (5.10) 
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Definition 5.1. We say that the perturbed problem has property EE(k) for some k > 3 if 
there exist S, t , C > and a function F : [0, 5) — > M + satisfying F(z) < Cz for z G [0, 5) so 
that the following hold. For any fjo, Vq, do satisfying 

\\(Vo,v ,a )\\ Hk < 5 (5.11) 

there exist (rj,v,cr) G £°°((0, t ); if 3 (f2)) so that 

1. (fj,v,a)(0) = (f) ,v ,o- ), 

2. r](t) = Id + fj(t) is invertible and 77 _1 (t) = Id — ((t) for < t < to, 

3. fj,v,a,( solve the perturbed problem on Q x (0,to), and 
4- we have the estimate 

sup \\(fj,v,a)(t)\\ H3 <F(\\(fj ,v ,a )\\ Hk ). (5.12) 
o<t<t 

Here the EE stands for existence and estimates, i.e. local-in-time existence of solutions 
for small initial data, coupled to L°°((0, to); H 3 (Q)) estimates in terms of the H k (Q) norm 
of the initial data. If we were to add the additional condition that such solutions be unique, 
then this trio could be considered a well-posedness theory for the perturbed problem. 

We can now show that property EE(k) cannot hold for any k > 3. The proof utilizes the 
Lipschitz structure of F to show that property EE(k) would give rise to certain estimates of 
solutions to the linearized equations ( 12. lip that cannot hold in general because of Theorem 

KB 

Theorem 5.2. The perturbed problem does not have property EE(k) for any k > 3. 

Proof. Suppose by way of contradiction that the perturbed problem has property EE(k) for 
some k > 3. Let 5, to, C > and F : [0, 5) — > 1R + be the constants and function provided by 
property EE(k). Fix n e N so that n > C. Applying Theorem 14.51 with this n, To = to/2, 
k > 3, and a = 1, we can find f/, v, a solving (12. lip so that 

\\^v,aW)\\m<- (5-13) 

but 

11^)11^3 > 11^)11^3 >1 for t> t /2. (5.14) 

For e > we then define fj^ = efj(0), = ev(0), and <7q = ea(0). 

Then for e < 5n we have || (fj^, v^, cXoOllfffc < ^> so according to EE{k) there exist 
(f/ e , v £ , a e ) G L°°((0, to); H 3 (fl)) that solve the perturbed problem with (^o>^o;^o) as m *~ 
tial data and that satisfy the inequality 

sup m,v^a £ )(t)\\ H3 < F(\\(f ,v s ,<j*)\\ Hh ) < Ce IKMaXO)^ < e. (5.15) 
o<t<t 
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Now define the rescaled functions if = r] £ /e,v £ = v £ /e,a £ = <j £ /e; rescaling (I5.15P then 
shows that 

sup \\(f,v £ ,a £ )(t)\\ H3 <1. (5.16) 
o<t<t 

Note that by construction (fj e , v £ , o" e )(0) = (77, v, er)(0). Our goal is to show that the rescaled 
functions converge as e — > to the solutions (fj, v, a) of the linearized equations (12. lip . 
We may further assume that e is sufficiently small so that 

sup \\ea £ {t)\\ Loo < - inf p (i 3 ) (5.17) 

0<t<t * -m<x 3 <£ 

and e < 1/(2K\), where K\ > is the best constant in the inequality 

||FG||„ 2 < i^FH^ ||G||„ 2 (5-18) 

for 3x3 matrix- valued functions F, G. The former condition implies that po + ea £ is bounded 
above and below by positive quantities, whereas the latter guarantees that 

B £ ■= (I- (I + eDfj T )- 1 )/e (5.19) 

is well-defined and uniformly bounded in L°°((0,to); H 2 (Q)) since 



B 



e 



H2 n=l 



E(-^rW)" 

n=l 

< j^r 1 wmi* < E < E = 2 - ( 5 - 2 °) 

n=l n=l n=l 

Since Id + efj £ is invertible, we may define ( £ via (Irf + e^) -1 = Id — e( £ , which implies that 
= r7 e o (Id - £C £ ). The slip map 5f. : JR 2 x K+ R 2 x {0} is then given by 

Si = Id R2 + ef + - e£ o (Jd Ra + e^). (5.21) 

The bounds on fj £ and the equation satisfied by ( £ then imply that 

sup \\Si(t) -Id R 2 II Loo <2e sup ||r7 £ (t)|| Loo < 2eAT 2 sup \\fj E (t)\\ H3 < 2eK 2 , (5.22) 

0<t<t o 0<<<t 0<t<t 

where K 2 > is the embedding constant for the trace map H 3 (Q) <- » L°°(R 2 x {0}). This 
bound allows us to define the normalized slip map Sf. := (Si —Idyp)/e as a well-defined and 
uniformly bounded function in L°°((0, to); L°°(M. 2 x {0}). Finally, we define the normalized 
remainder function by 

rEo- E (x,t) 



TZ £ (x, t) = - (ea £ (x, t) - z)h"(p (x) + z)dz 

£ Jo 



r<r E (x,t) 

el (a £ (x,t) - z)h"(p (x) +ez)dz. (5.23) 
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A straightforward calculation employing (I5.16P and (I5.17P shows that 

sup ||^ £ 0)|L 3 < sK 3 (5.24) 

0<t<to 

for some constant K 3 > 0. 

We can now parlay the bounds on f] £ , v £ , a £ , B £ , and TZ £ into corresponding bounds on 
dtf] e , dtv £ , and dta £ and some convergence results. The first equation in (15. 5p implies that 
d t fj e = v £ , so that 

sup \\d t fj £ {t)\\ H3 = sup 11^)11^3 <1. (5.25) 

0<t<t 0<t<t 

Expanding the second equation in (15. 5p . we find that 

(d t a £ + po div v £ ) + e(a £ div v £ - p tr(B £ Dv £ )) - s 2 (a £ tr(B £ Dv £ )) = 0. (5.26) 
This implies that 

lim sup \\d t a £ (t) + p divv £ (t)\\ H2 =0 (5.27) 

£ ^°0<t<t 

and that 

sup \\d t cr £ (t)\\ H2 < K 4 (5.28) 

o<t<t 

for a constant K4 > 0. Performing a similar expansion on the third equation in (15. 5p yields 

(d t v £ + V(h'( Po )a £ + ge 3 ■ rf)) + e(VTZ £ - B £ V(h'(p )a £ + ge 3 ■f) £ )) + e 2 (B £ VTZ £ ) = 0. (5.29) 
Hence 

lim sup \\d t v £ (t) + V(h'(p )a £ (t)+ge 3 -r ] £ (t))\\ H2 =0 (5.30) 

o<t<t 

and 

sup \\d t v £ (t)\\ H2 <K 5 (5.31) 
o<t<t 

for a constant K§ > 0. 

We now turn to some convergence results for the jump conditions. We begin with the 
second equation in (15.71) . which we expand using the mean- value theorem to get 

P + (p+) + P|«Po + + (1 - of+)ea + )ed% 

= P_(p ) + P'_(a £ _ Po + (1 - al_)ea„ o (Jd R3 + eS £ _))eat o (Jd Ra + eS £ _) (5.32) 

for functions a £ ± : R 2 x IR + — > [0, 1]. From the above bounds on a £ and S £ we have that 

sup ||P;«p + + (1 - a £ + )ea + ) - PM)\\ LOO -> 0, (5.33) 

0<t<t o 

sup ||p:(«Ip + (1 - (^k 2 + ^)) " ^-(Po )|U 0, (5.34) 

o<t<t 

and 

sup IIctI o (Jd R 2 + eS"!) - alll < sup || Va £ (t)\\ Lao sup lleS^) II as e 0. 

o<t<t o<t<t o<t<t 

(5.35) 
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Since P+(Pq ) = P-(Po )> we ma y eliminate these terms from equation (15.321) and divide both 
sides by e; then employing (I5.33I) - (I5.35I) . we deduce that 

sup ||P|(p+)a;(t) - PL{Po)dt{t)\\ Laa -> as e 0. (5.36) 

0<t<t 

For the second equation in (15 .7p we first write the normal at the interface as n £ = N £ / \N £ \ 
with 

N £ = (ex + ed Xl f] e + ) x (e 2 + ed X2 fj E + ) 

= e 3 + e(ei x 9 xa ^. + ^77* x e 2 ) + £ 2 (<9 xi r^ x d^) := e 3 + eiV £ . (5.37) 

As e — > we have |iV e | > 0, so we may rewrite the second equation in (15.71) as 

(v £ + - vl o (Id R 2 + eS e _)) ■ (e 3 + eN £ ) = 0. (5.38) 

Clearly sup 0<t<<0 ||^ e (^)|| i0 o is bounded uniformly and 

sup \\v £ _o(Idu2+eS £ _)—v l L\\ < sup ||D-u e (t)|| LOO sup | \eS e (t) | — > as e — >■ 0. 

o<t<t o<t<t o<t<t 

(5.39) 

Using these, we find that 

sup ||e 3 • (v £ + (t) -v £ _(t))\\ as e 0. (5.40) 

o<t<t 

According to the bound (15.161) and sequential weak-* compactness, we have that up to 
the extraction of a subsequence (which we still denote using only e) 

[f] £ ,v £ ,a £ ) ^ (f)°,v°,a ) weakly-* in L°°((0, t ); H 3 {Q)). (5.41) 

By lower semicontinuity we know that 

sup \\(f} ,v ,a°)(t)\\ H3 <l. (5.42) 

0<t<t 

On the other hand, by (I5T25D . fl53H|) . and floTBTD . we know that 

limsup sup \\(d t f] £ ,dtV £ ,d t a £ )(t)\\ H 2 < oo. (5.43) 

£->0 o<t<t 

By a result in [6], we then have that the set {(f] £ , v £ , a £ )} is strongly pre-compact in the 
space L°°((0,to);# 11/4 (fi)), so 

(rf,v £ ,a £ ) (r/°,w°,a ) strongly in L°°((0, t ); H ll/ \ti)). (5.44) 

This strong convergence, together with the convergence results (I5.27p . (I5.30P and the equa- 
tion d t ff = v £ , implies that 

(d t fj £ ,d t v £ ,d t a £ ) -> (d t f} ,d t v°,d t a ) strongly in L°°((0, t ); # 7/4 (^)), (5.45) 

29 



and that 

d t ff = v° 

d t a° + podivv = (5.46) 
d t v° + V{h'{p )a° + ge 3 ■ ff) = 0. 

We may pass to the limit in the initial conditions (fj e , v 6 , <x £ )(0) = (fj, v, fr)(0) to find that 

(fj°,v°,a )(0) = (fj,v,a)(0) (5.47) 

as well. We now derive the jump and boundary conditions for the limiting functions. The 
index 11/4 is sufficiently large to give L°°((0,t ); L°°) convergence of (f] e ,v £ ,a £ ) when re- 
stricted to {x 3 = 0}, {x 3 = —m}, and {x 3 = £}, i.e. the interface and the lower and upper 
boundaries. Combining this with (15.361) and (15.401) . we deduce that 

P'M)°l = P-(Po)°°- on {x 3 = 0}, (5.48) 

(v° + -v°_) -e 3 = on {x 3 = 0}, (5.49) 

and that 

v+ ■ e 3 = on {x 3 = £} and v°_ ■ e 3 = on {x 3 = —m}. (5.50) 

Now, according to (I5.46I) - (I5.50I) . (fj , v°, a ) are solutions to 02. 1 lj) and the corresponding 
jump and boundary conditions on Q x (0, to) that satisfy the initial condition (I5.47p . Then 
according to Theorem I4.4[ 

(fj°, v°, a ) = (fj, v, a) on x [0, t ). (5.51) 
Hence we may chain together inequalities (I5.42p and (I5.14p to get 

2< sup IKr/W^WLa < sup ||(77 ,^ ,o-°)(t)||^3 < 1, (5.52) 

t /2<t<t ' 0<t<t 

which is a contradiction. Therefore, the perturbed problem does not have property EE(k) 
for any k > 3. □ 
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